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WAVEFRONT CODING INTERFEFENCE CONTRAST IMAGING SYSTEMS 
BACKGROUND OF THE INVENTION 
FIELD OF THE INVENTION : 

This invention relates to apparatus and metliods for using Wavefront Coding 
to Improve contrast Imaging of objects wiiich are transparent, reflective or 
vary in tliiclcness or index of refractiori. 

DESCRIPTION OF THE PRIOR ART: 

U.S. Patent No. 5,748,371, issued May 5, 1998 and entitled "Extended Deptli 
of Field Optical Systems," is a relevant reference. 

Most imaging systems generate Image contrast through variations in 
reflectance or absorption of the object being viewed. Objects that are 
transparent or reflective but have variations in index of refraction or 
thicl<ness can be very difficult to image. These types of transparent or 
reflective objects can be considered "Phase Objects". Various techniques 
. have been developed to produce high contrast images from essentially 
transparent objects that have only variations in thickness or index of 
refraction. These techniques generally modify both the illumination optics 
and the imaging optics and are different modes of what can be called 
"Contrast Imaging". 

There are a number of different Contrast Imaging techniques that have been 
developed over the years to image Phase Objects. These techniques can be 
grouped into three classes that are dependent on the type of modification 
made to the bacl< focal plane of the imaging objective and the type of 
illumination method used. The simplest Contrast Imaging techniques modify 
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the back focal plane of the imaging objective with an Intensity or amplitude 
mask. Other techniques modify the back focal plane of the objective with 
phase masks. Still more techniques require the use of polarized illumination 
and polarization-sensitive beam splitters and shearing devices. 

5 Contrast Imaging techniques that require polarizers, beam splitters and 

beam shearing to image optical phase gradients, we call "Interference 
Contrasf techniques. These techniques include conventional Differential 
Interference Contrast (Smith, L.W., Microscopic Interferometry, Research 
(London), 8:385-395, 1955), improvements using Nomarski prisms (Allen, R. 

I 0 D., David, Q. B, and Nomarski, G, The Zelss-Nomarski differential interference 

equipment for transmitted light microscopy, Z. Wiss. Mikrosk. 69:193-221, 
1969), the Dyson interference microscope (Bom and Wolf, Principals of 
Optics, Macmlllan, 1964), the Jamin-Lebedeff Interferometer microscopes as 
described by Spencer in 1982 ("Fundamentals of Light Microscopy", 

1 5 Cambridge University Press, London), and Mach-Zehnder type Interference 

microscopes ("Video Microscopy", Inoue and Spring, Plenum Press, NY, 
1997). Other related techniques include those that use reduced cost beam 
splitters and polarizers (US Patent # 4,964,707), systems that employ 
contrast enhancement of the detected images (US Patent # 5,572,359), 

2 0 systems that vary the microscope phase settings and combine a multiplicity 

of images (US Patent # 5,969,855), and systems having variable amounts 
of beam shearing (US Patent # 6,128,127). 

Figure 1 (Prior Art) Is a block diagram 100, which shows generally how 
Interference Contrast Imaging techniques are implemented. This block 
25 diagram shows Imaging of a Phase Object 110 through transmission, but 

those skilled in the art will appreciate that the elements could just as simply 
have been arranged to show imaging through reflection. 

Illumination source 102 and polarizer 104 act to fom) linearly polarized light. 
Beam splitter 106 divides the lineariy polarized light Into two llneariy 
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polarized beams that are orthogonally polarized. Such orthogonal beams 
can be laterally displaced or sheared relative to each other. Illumination 
optics 108 act to produce focussed light upon Phase Object 110. A Phase 
Object is defined here as an object that is transparent or reflective but has 
5 variations in thickness and/or Index of refraction, and thus can be difficult to 

image because the majority of the image contrast typically is derived from 
variations in the reflectance or absorption of the object. 

Objective lens 112 and tube lens 118 act to produce an image upon 
detector 120. Beam splitter 114 acts to remove the lateral shear between 

1 0 the two orthogonally polarized beams formed by beam splitter 106. Beam 

splitter 114 is also generally adjustable. By adjusting this beam splitter a 
phase difference between the two orthogonal beams can be realized. 
Analyzer 116 acts to combine the orthogonal beams by converting them to 
the same linear polarization. Detector 122 can be film, a CCD detector 

1 5 array, a CMOS detector, etc. Traditional imaging, such as bright field 

imaging, would result if polarizer and analyzer 104 and 116 and beam 
splitters 106 and 114 were not used. 

Figure 2 (Prior Art) shows a description of the ray path and polarizations 
through the length of the Interference Contrast imaging system of Figure 1. 

20 The lower diagram of Figure 2 describes the ray path while the upper 

diagram describes the polarizations. The illumination light is linearly 
polarized after polarizer 204. This linear polarization is described as a 
vertical arrow in the upper diagram directly above polarizer 204. At beam 
splitter 206 the single beam of light becomes two orthogonally polarized 

25 beams of light that are spatially displaced or sheared with respect to each 

other. This is indicated by the two paths (solid and dotted) in both 
diagrams. Notice that the two polarization states of the two paths in the 
top diagram are orthogonally rotated with respect to each other. Beam 
splitter 214 spatially combines the two polarizations with a possible phase 
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offset or bias. This phase bias is given by the parameter A in the upper 
plot. By laterally adjusting the second beam splitter 214 the value of the 
phase bias A can be changed. A Nomarski type prism is described by the 
ray path diagram, although a Wollaston type prism could have been used as 
well. Analyzer 216 acts to convert the orthogonal component beams to 
linearty polarized light. The angle between the polarizer 204 and analyzer 
216 can typically be varied in order to adjust the background intensity. 
Image plane 21 8 acts to display or record a time average intensity of the 
linearly polarized light, the sheared component possibly containing a phase 
shift. This image plane can be an optical viewing device or a digital detector 
such as CCD, CMOS, etc. 

The interactions of the polarizers, beam splitters, and Phase Objects of the 
Interference Contrast imaging systems have been studied in great detail. For 
additional background information see "Confocal differentia! interference 

1 5 contrast (DIG) microscopy: including a theoretical analysis of conventional 

and confocal DIC imaging", Cogswell and Sheppard, Journal of Microscopy, 
Vol 165, Pt 1, January 1992, pp 81-101. 

In order to understand the relationship between the object, image, and 
phase shift A consider an ari^itrary spatially constant object that can be 

2 0 mathematically described as: 

Obj ^ a exp( je), where j = V- 1 

where "a" is the amplitude and e is the object phase. If the two component 
beams of the system of Figure 2 have equal amplitude, and if the 
component beams are subtracted with relative phases +/- A/2 then just 
2 5 after analyzer 216 the resulting image amplitude is given by: 

amp = a exp( ][e-A/2] ) - a exp( j[e+A/2] ) 2 j a exp( je ) sin(A/2) 

The image intensity is the square of the image amplitude. The Intensity of 

4 
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this signal is tlien given by: 

int^ = 4 a2 sin(A/2)2 

The image intensity is independent of the object phase e. The phase 
difference or bias between the two orthogonal beams is given by A and is 
5 adjusted by lateral movement of the beam splitter, be it a Wollaston or a 

Nomarski type. If instead of a spatially constant object, consider an object 
whose phase varies by A* between two laterally sheared beams. This 
object phase variation is equivalent to a change in the value of the 
component beam phases of A. If the component beam phases A is equal to 

1 0 zero (no relative phase shift) then the resulting image intensity can be 

shown to have Increases in intensity for both positive and negative 
variations of object phase. If the component beam bias is increased so 
that the total phase variation Is always positive, the change in image 
intensity then increases monotonically throughout the range AO. The actual 
15 value of the change in image intensity with object phase change AO can be 

shown to be: 

lnti = 4 a2 AO sin(A) 

In Interference Contrast imaging the phase bias A determines the relative 
strengths with which the phase and amplitude infomnation of the object will 

2 0 be displayed in the image. If the object has amplitude variations these will 

be imaged according to int^ above. ^ At a phase bias of zero (or multiple of 
2pi ) the image will contain a maximum of phase infomnation but a minimum 
of amplitude information. At a phase bias of pi the opposite is true, with 
the image giving a maximum of amplitude information of the object and a 
25 minimum of phase information. For intermediate values of phase bias both 

phase and amplitude are imaged and the typical Interference Contmst bas 
relief image is produced, as is well known. 

5 
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Variation of the phase bias can be shown to effect the parameters of 
image contrast, linearity, and signal-to-noise ratio (SNR) as well. The ratio 
of contrast from phase and amplitude in Interference Contrast imaging can 
be shown to be given by: 

[ contrast due to phase / contrast due to amplitude ] = 2 cot(A/2) 

The overall contrast In the Interference Contrast image is the ratio of the 
signal strength to the bacl<ground and can be shown to be given by: 

overall contrast ^ 2 AO cot(A/2) 

The linearity between the image intensity and phase gradients in the object 
can be described by: 

L = [ ( 1+sin(A))(2/3) ] / [ 2 cos(A) 1 

The signal-to-noise ratio (SNR), ignoring all sources of noise except shot 
noise on the background, can be shown to be given by 

SNR = 4 a cos( A/2 ) 

In Interference Contrast imaging systems the condenser aperture can be 
opened to improve resolution, although in practice, to maintain contrast, the 
condenser aperture is usually not increased to full illumination. Imaging is 
typically then partially coherent. Description of the imaging characteristics 
for Interference Contrast imaging therefore needs to be expressed in terms 
of a partially coherent transfer function. The partially coherent transfer 
function (or transmission cross-coefficient), given as C(m,n;p,q), describes 
the strength of image contributions from pairs of spatial frequencies 
components m; p in the x direction and n; q in the y direction (Bom and 
Wolf, Principals of Optics, Macmillan, 1975, p. 526). The intensity of the 
image in terms of the partially coherent transfer function image can be 
written as: 
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l(x,y) = J 1 n T(m,n) T(p,g)*C(m,n;p,q) 

exp(2 pi j [(m-p)x+(n-q)y] ) dm dn dp dq 

where the limits of Integration are + infinity to - infinity. The term T(m,n) Is 
the spatial frequency content of the object amplitude transmittance t(x,y): 

T(m,n) = J J t(x,y) exp( 2 pi j [mx+ny] ) dx dy 

where again the limits of integration are + infinity to - infinity, ( )* 
denotes complex conjugate. When the condenser aperture is maximally 
opened and matched to the bac!< aperture or exit pupil of the objective lens, 
the partially coherent transfer function reduces to (Intro, to Fourier Optics, 
Goodman ,1968, pg.120): 

C(m, n; p, q) = 6(m-n) 5(p-q) [ acos(p) - p sqrt{(1-p2)} ] 

where p = sqrt( m2+ p2 ) and 8{x) = T if x=0, 5(x) = 0 otherwise. 

The effective transfer function for the Interference Contrast imaging system 
can be shown to be given as: 

C{m,n;p,q)3ff = 2 C(m,n;p,q) { cos[2 pi(m-n)A] - 

cos(A) cos([2 pi(m+n) A] - sin(A)sin[2 pi (m+p) A] } 

where A is equal to the lateral shear of the beam splitters and C(m,n;p,q) is 
the partially coherent transfer function of the system without Interference 
Contrast modifications. 

Interference Contrast imaging is one of the most complex fomns of imaging 
in terms of analysis and design. These systems are also widely used and 
studied. But, there is still a need to improve Interference Contrast Imaging 
of Phase Objects by increasing the depth of field for imaging thicl< objects, 
as well as for controlling focus-related aberrations in order to produce less 
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expensive imaging systems than is currently possible. 

SUMMARY OF THE INVENnON 

An object of the present Invention is to improve Contrast Imaging of Phase 
Objects by increasing depth of field and controlling focus-related 
aberrations. This is accomplished by using Contrast Imaging apparatus and 
methods with Wavefront Coding aspheric optics and post processing to 
increase depth of field and reduce misfocus effects. The general 
Interference Contrast imaging system is modified with a special purpose 
optical element and image processing of the detected image to form the 
final image. Unlike the conventional Interference Contrast Imaging system, 
the final Wavefront Coding Interference Contrast image Is not directly 
available at the Image plane. Post processing of the detected image is 
required. The Wavefront Coding optical element can be fabricated as a 
separate component, can be constmcted as an integral component of the 
Imaging objective, tube lens, beam splitter, polarizer or any combination of 
such. 

Apparatus for increasing depth of field and controlling focus related 
aberations in an Interference Contrast Imaging system having an 
illumination source, optical elements for splitting light polarizations, and 
20 illumination optics placed before a Phase Object to be imaged, and elerrlents 

for recombining light polarizations and objective optics after the Phase 
Object to form an image at a detector, includes an optical Wavefront 
Coding element having an aperture and placed between the Phase Object 
and the detector, the coding element being constructed and arranged to 
25 alter the optical transfer function of the Interference Contrast Imaging 

system in such a way that the altered optical transfer function is 
substantially insensitive to the distance between the Phase Object and the 
objective optics over a greater range of object distances than was 
provided by the unaltered optical transfer function, wherein the coding 

8 
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element affects the alteration to the optical transfer function substantially 
by affecting the phase of light transmitted by the element. The system 
further includes a post processing element for processing the image 
captured by the detector by reversing the alteration of the optical transfer 
function accomplished by the coding element. 

The detector might be a charge coupled device (CCD). 

The phase of light transmitted by the coding element is preferably relatively 
flat near the center of the aperture with increasing and decreasing phase 
near respective ends of the aperture. 

As an alternative, the phase of light transmitted by the coding element could 
substantially follow a cubic function. 

In one embodiment, the phase of light transmitted by the coding element 
substantially follows a function of the form: 

Phase (x,y) = 12[x^+y^] 

where Ixl ^ 1 , lyl ^ 1 . 

In another embodiment the phase of light transmitted by the coding element 
substantially follows a rectangularly separable sum of powers function of 
the form: 

phase(x,y) = [ S a, sign{x) jxl C/Sign(y) lyl ^/] 

where the sum is over the index /, 

sign(x) = -1 for x < 0, sign(x)=+1 for x > 0. 

In another embodiment, the phase of light transmitted by the coding 
element substantially follows a non-separable function of the form: 
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phase(r,e) = [2 r^; cos( 6 + ) ] 
where the sum is again over the index /. 

In another embodiment the phase of light transmitted by the coding element 
substantially follows a function of the fomn: 

Phase profile (x,y) = 

7[ sign(x) lx|3 + sign(y) lyP ] + 7[ sign(x) Ixl^ ^ + sign{y) lyl^-^] 

where Ixi ^1, lyl ^ 1 

The coding element further may be integrally fonned with a lens element for 
focussing the light, or with the illumination optics. 

The coding element could comprise an optical material having varying 
thickness, an optical material having varying index of refraction, spatial light 
modulators, or micro-mechanical mirrors. 

A method for increasing depth of field and controlling focus related 
aberrations in a conventional Interference Contrast Imaging system 
comprises the steps of modifying the wavefront of transmitted light 
between the Phase Object and the detector, the wavefront modification 
step selected to alter the optical transfer function of the Interference 
Contrast Imaging system in such a way that the altered optical transfer 
function is substantially insensitive to the distance between the Phase Object 
and the objective optics over a greater range of object distances than was 
provided by the unaltered optical transfer function, and post processing the 
image captured by the detector by reversing the alteration of the optical 
transfer function accomplished by the element. 

A Wavefront Coding optical element can also be used on the Illumination side 
of the system in order to extended the depth of field of the projected 
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illumination due to the duality of projection and imaging. Tliis projected 
illumination would be broader than without Wavefront Coding, but the optical 
density as a function of distance from the object would be less sensitive 
with Wavefront Coding than without. Without Wavefront Coding on the 
Illumination side of the system, the object can technically be Imaged clearly 
but Is not illuminated sufficiently. See "Principal of Equivalence between 
Scanning and Conventional Optical Imaging Systems", Dorian Kermisch, J. 
Opt; Soc. Am., Vol. 67, no. 10, pp.1 357-1 360 (1977). 

BRIEF DESCBIPTIQN OF THE DRAWINGS 

Figure 1 (prior art) shows a standard prior art Interference Contrast 
imaging system. 

Figure 2 (prior art) shows ray paths and polarization states for the 
Interference Contrast imaging system of Figure 1 . 

Figure 3 shows a Wavefront Coding Interference Contrast imaging system 
including Wavefront Coding optics and post processing in accordance with 
the present invention. 

Figure 4 describes in detail the Object IVIodifylng Function and Object Imaging 
Function of the Wavefront Coding Interference Contrast system. 

Figure 5 shows the aperture transmittance function and the corresponding 
ambiguity function for the Object Imaging Function of the prior art system 
of Figure 1. 

Figure 6 shows the Wavefront Coded cubic phase function and the 
corresponding ambiguity function for the Object Imaging Function of Figure 
3. 

Figure 7 shows another Wavefront Coded phase function and the 
conesponding ambiguity function for the Object Imaging Function of Figure 
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3. 

Figure 8 shows misfocus MTFs for the prior art Object imaging Function of 
Figure 1 and the Object Imaging Functions for the Wavefront Coded 
interference Contrast systems described in Figures 3, 6 and 7. 

Figure 9 shows singie plane of focus images of human cervical cells with 
darkly stained nuclei imaged with a 40X. NA=1.3 objective with a 
conventional Interference Contrast system and with a Wavefront Coded 
interference Contrast imaging system similar to that of Figure 3. 

nPTAILED DESCP IPnONl QFTHE PRFFFRRFD EMBOPIMENTTS 

Wavefront Coding can be used with conventional objectives, polarizers and 
beam splitters In Interference Contrast systems, as shown in Figure 3, to 
achieve an increased depth of field in an optical and digital imaging system. 
This can be explained by considering the Object Modifying Functions of 
conventional Interference Contrast systems separately from the Object 
Imaging Functions, as shown In Figure 4. By considering these two functions 
separately, modification of depth of field can be explained in terms of the 
Object imaging Function. Extending the depth of field of the Object Imaging 
Functions of Interference Contrast systems is shown in Figures 5-8, Figure 
9 shows real-world images of human cen/icai cells taken with a system 
having only Interference Contrast and a comparison to an image from a 
Wavefront Coding Interference Contrast system. 

Figure 3 shows a Wavefront Coded Interference Contrast imaging system 
300 including Wavefront Coding and post processing in accordance with the 
present invention. Similar reference numbers are used in Figure 3 as are 
used in Figure 1 , since the systems are very similar, except for the addition 
of Wavefront Coding element 324 and post processing 326. The general 
Interference Contrast imaging system of Figure 1 is modified with a special 
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purpose generalized aspheric optical element 324 and image processing 326 
of the detected image to form the final image, Unlli<e the conventional 
Interference Contrast system, the final image In combined system 300 is not 
directly available at detector 322. In fact, no sharp and clear image of any 
kind is available in system 300, except after image processing 326. Image 
processing 326 of the detected image Is required to remove the spatial 
Wavefront Coding effects (other than the extended depth of field). 

Wavefront Coding optical element 324 can be fabricated as a separate 
component as shown In figure 3, or can be combined with objective lens 
312, tube lens 318, beam splitter 314, analyzer 316, or any combination of 
these. Any material or configuration that can impart a range of spatial 
phase shifts to a wavefront can be used to construct Wavefront Coding 
element 324. For example, optical glass or plastic of varying thickness 
and/or index of refraction can be used. Holograms and mirrors can also be 
used as the material for the Wavefront Coding element. In order to 
dynamically adjust the amount of depth of field, or to eissentially change the 
Wavefront Coding element 324 for different objectives or desired depth of 
field, spatial light modulators or dynamically adjustable micro mirrors or 
similar can also be used. 

Wavefront Coding optical element 324 can also be used on the illumination 
side of system 300 In order to extend the depth of field of the projected 
Illumination due to the duality of projection and imaging. This projected 
illumination would be broader than without Wavefront Coding, but the optical 
density as a function of distance from the object would be less sensitive 
with Wavefront Coding than without. 

The components that distinguish the Wavefront Coding Interference 
Contrast system of Figure 3 from a general or brightfield imaging system is 
polarizer 304, beam splitter 306, beam splitter 314, analyzer 316, and 
Wavefront Coding element 324 and Image processing 326. The polarizer, 
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analyzer, and beam splitters essentially use phase to modify the imaging 
characteristics of the object 310. The Wavefront Coding element 324 and 
image processing 326 are used to increase the depth of field or remove 
misfocus effects in images of the modified object as shown below. By 
grouping the components of system 300 by their function, the Wavefront 
Coding Interference Contrast Imaging system can be understood. 

The locations of polarizer, analyzer, and beam splitters of Figure 3 have 
been chosen because of historical reasons. These are tlie traditional 
locations for these components in prior art systems relative to the 
illumination and imaging optics. The same relative locations are seen in 
Figure 1. The beam splitter 314 and analyzer 316 can theoretically be 
moved relative to objective lens 312 without changing the imaging behavior 
of the system. See system 400A of Figure 4. Numbering conventions of 
Figure 4 are also similar to those of Figures 1 and 3 due to the similar 
nature of the components, in system 400A the beam splitter and analyzer 
have been moved before the objective lens but after the object The 
wavefront after analyzer 416 is linearly polarized as is the wavefront after 
analyzers 216 and 316 in Figures 2 and 3 respectively. Since, ideally, lenses 
do not change the polarization, shear, or bias of the wavefront this new 
location is technically equivalent to that of Figure 3. Consider the ray paths 
of Figure 2. Notice that the ray paths between beam splitters 206 and 214 
are parallel. Moving beam splitter 206 before objective lens 212 
theoretically will not ctiange the parallel nature of the ray paths. Analyzer 
216 can also move before objective lens 212 with no adverse affects. The 
component arrangement of system 40aA allows the "Object Modifying 
Functions" to be cleariy distinguished from the Object Imaging Functions. 

in order to further characterize the Object Modifying Function of system 
400A consider system 400B of Figure 4. In this system a new phase and 
amplitude object 4108 replaces the original object 41 OA of system 400A. 
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This new object is selected so tiiat its three dimensional structure produces 
an identical wavefront from Illumination source 402, polarizer 404, and 
Illumination optics 408 as from object 41 OA when combined with the 
polarizer, analyzer, and beam splitters of system 400A. It is well known 
that a phase and amplitude object can be theoretically constructed so that 
any given linearly polarized wavefront can be reproduced from linearly 
polarized illumination. Although it is theoretically possible to produce such a 
new object 41 OB, in practice it might be difficult. Since a new object 41 OB 
can be substituted for the combination of original object 41 OA, beam 
splitter 406, beam splitter 414, and analyzer 416, it is clear that the 
polarizers and analyzers do nothing more than simply act to modify the 
imaging characteristics of the object. Notice that the right sides of 
systems 400A and 400B are identical. The right sides of these systems are 
the Object Imaging Function. The Object Imaging Function images the 
object that has had its imaging characteristics modified by the Object 
l\/lodifying Function. With the Wavefront Coding optical element 424 and 
image processing 426 the Object Imaging Function can have a very large 
depth of field and be able to control focus-related aberrations. 

If the Object Imaging Function of system 400B has a large depth of field, 
then the New Object of 41 OB can be imaged over a large depth. Likewise, 
when the Object imaging Function of figure 400A has a large depth of field, 
object 41 OA (as modified by the Object Modifying Function) can be imaged 
with a large depth of field. Since system 400B produces identical images to 
system 400A, and system 400A produces identical images to system 300, 
this also means that system 300 will image object 310 with a large depth of 
field. This large depth of field is also independent of the object or Object 
Modifying Functions as shown in Figure 4. 

The Object Imaging Function can be made to have a large depth of field by 
use of a generalized aspheric optical element and signal processing of the 
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detected images. Ambiguity function representations can be used to 
succinctly describe tinis large depth of field. Only the magnitude of the 
ambiguity functions in this and following figures are shown. Ambiguity 
functions are, In general, complex functions. One-dimensional systems are 
given for simplicity. Those skilled in the art of linear systems and ambiguity 
function analysis can quickly make extensions to two-dimensional systems. 
An ambiguity function representation of the optical system is a powerful 
tool that allows MTFs to be Inspected for all values of misfocus at the same 
time. Essentially, the ambiguity function representation of a given optical 
system is similar to a polar plot of the MTF as a function of misfocus. The 
in-focus MTF is described by the trace along the horizontal v=0 axis of the 
ambiguity function. An MTF with normalized misfocus value of y = 2 pi X. 
W20, where W20 is the traditional misfocus aberration coefficient and X is the 
illumination center wavelength, is described in the ambiguity function along 
the radial line with slope equal to (\|f / pi). For more infonmation on 
ambiguity function properties and their use in Wavefront Coding see 
"Extended Depth of Field Through Wavefront Coding", E. R. Dowski and W.T. 
Cathey, Applied Optics, vol. 34, no 11, pp. 1859-1866, April. 1995. and 
references contained therein. 

Figure 5 gives an ambiguity function perspective on the Object Imaging 
Function of conventional Interference Contrast systems. The top plot of 
Figure 5 shows the aperture transmittance function of an ideal conventional 
Interference Contrast system such as that shown in Figure 1. The bottom 
plot of Figure 5 shows the associated ambiguity function associated with 
the Object Imaging Function for the prior art system of Figure 1 . 

Over the nonnalized aperture (in nomiallzed coordinates extending from -1 
to +1) the conventional system has a transmittance of 1, i.e. 100%, The 
phase variation (not shown) is equal to zero over this range. The 
corresponding ambiguity function has concentrations of optical power 
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(shown as dark shades) very close to the horizontal v=0 axis. From the 
relationship between the ambiguity function and misfocused MTFs, we see 
that the conventional interference Contrast Systems has a small depth of 
field because slight changes In misfocus lead to MTFs (represented by radial 
5 lines with non-zero slope in the ambiguity function) that intersect regions of 

small power. 

Figure 6 shows an example of a phase function for the Wavefront Coding 
optical element 324 and corresponding ambiguity function for an improved 
system of Figure 3. This phase function is rectangularly separable and can 
10 be mathematically described in two dimensions as: 

Phase (x,y) = 12 [x^+yS] 
ix I ^ 1 , lyl ^ 1 

Only one dimension of this phase function Is shown in the upper plot of Figure 
6. Increasing the peak-to-valley phase height (as can be done by increasing 
1 5 the constant 12 above) results in increasing depth of field. The 

transmittance of this system (not shown) is unity (I.e. 100%) over the entire 
aperture, as in the top plot of Figure 5. 

The ambiguity function shown In Figure 6 for this Wavefront Coded 
Interference Contrast system is seen to have optical power spread over a 
much larger region in the ambiguity domain than does the diffraction-limited 
system plotted in Figure 5. Broader regions of optical power in the 
ambiguity function translate to larger depth of field or depth of focus since 
the ambiguity function is essentially a radial plot of misfocused MTFs with 
the angular dimension pertaining to misfocus. Thus, this Wavefront Coded 
Interference Contrast system has a larger depth of field than the 
conventional interference Contrast system. 

There are an infinite number of different Wavefront Coding phase functions 
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that can be used to extend the depth of field. Other more general 
rectangularly separable forms of the Wavefront Coding phase function is 
given by: 

phase{x,y) = [ S a/Sign(x) Ixl ^+ CySign(y) lyl ^j] 
where the sum is over the index /, 
slgn(xj = -1 for x < 0, slgn(x)=+1 for x ^ 0. 

Rectangularly separable forms of Wavefront Coding allow fast processing. 
Other ioms of Wavefront Coding complex phases are non-separable, and 
the sum of rectangularly separable forms. One non-separable form is 
defined as: 

phase{r,e) = [Z r^i cos( by 6 + ) ] 

where the sum is again over the index /. 

Figure 7 shows the Wavefront Coding phase function and the ambiguity 
function for a further improved system of Figure 3. The top plot of Figure 7 
shows the phase function from Figure 6 (curve 701) and a further improved 
phase function (curve 702). The aperture transmittance function is the 
same as shown in Figure 5. The form of the new phase profile 702, in 
radians, of this system is given by: 

Phase profile (x,y) = 

7[ sign(x) lx|3 + sign(y) lyis ] + 71 sign(x) Ixl^^e + sign(y) lyl^ e ] 

1x1 ^1, lyl <. 1 

The ambiguity function related to phase function 702 is shown in the bottom 
of Figure 7. This ambiguity function is seen to have more optical power 
unifomily spread about the horizontal v=0 axis When compared to either the 
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Wavefront Coding Interference Contrast system plotted in Figure 6 or the 
Conventional Interference Contrast system plotted in Figure 5. Thus, the 
Wavefront Coded Interference Contrast system of Figure 7 will have a larger 
depth of field than the systems represented In Figures 6 or 5. 

Figure 8 shows modulation transfer functions as a function of misfocus for 
the prior art Interference Contrast system, and the Wavefront Coded 
Interference Contrast systems of Figures 6 and 7. The top plot of Figure 8 
shows the MTFs of the conventional Interference Contrast Imaging system 
of Figures 1 and Figure 5 and the MTFs of the Wavefront Coded Interference 
Contrast system of Figure 6. The bottom plot shows the MTFs of the 
Interference Contrast imaging system of Figures 1 and 5 (again) and the 
MTFs from the Wavefront Coding Interference Contrast imaging system of 
Figure 7. These plots are the particular MTFs given in the respective 
ambiguity functions for the normalized misfocus values \|f = { 0, 2, 4}. 
Notice that the MTFs for the conventional Interference Contrast system 
(top and bottom plots) vary appreciably with even this slight amount of 
misfocus. The image from the conventional system will thus change 
drastically due to misfocus effects for only small, misfocus values. This is 
expected from the narrow ambiguity function associated with the 
conventional system (shown In Figure 5). 

By comparison, the MTFs from the Wavefront Coded Interference Contrast 
imaging systems (top and bottom plots) show very little change with 
misfocus as predicted by the ambiguity functions associated with these 
systems (shown in Figures 6 and 7). If the MTFs of the system do not 
change, the resulting MTFs (and hence also point spread functions) can be 
corrected over a large range of misfocus with a single post processing step 
326. A single post processing step is not possible with conventional 
systems, which change appreciably with misfocus since the MTFs and PSFs 
of the system change with misfocus to values that are unknown and often 
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impossible in practice to calculate. The MTFs from the Wavefront Coded 
Interference Contrast system in the top plot are seen to have lower values 
for most spatial frequencies than the MTFs from the Wavefront Coded 
Interference Contrast system of the bottom plot. This is expected from the 
ambiguity functions of Figures 6 and 7 respectively. The two-temn phase 
function (curve 702) yields MTFs that not only have similarly small change 
with misfocus but also give a higher MTF then those associated with the 
simple cubic phase function (curve 701). This higher MTF results in a more 
compact PSF (not shown) as well as less signal-to-noise ratio penalties 
needed for the image processing 326. 

In general, the Wavefront Coded objective element phase function that yields 
the smallest MTF variation with misfocus and also the highest MTF is 
preferred in practice. There are an infinite number of different objective 
masl< phase functions that are good candidates for control of the MTF. The 
characteristics that practical Wavefront Coding element phase functions 
have can generally be described as being relatively flat near the center of 
the aperture with increasing and decreasing phase near the respective 
edges of the aperture. The central portion of the phase function controls 
the majority of the light rays that would not need modification If the 
objective were stopped down, for the depth of field extension required. For 
increasing amounts of depth of field, the size of the central phase region 
that can be fiat decreases. Increasing the flatness of the central region of 
the rays leads to larger MTFs as seen in comparison to the phase functions 
and MTFs of Figures 6, 7, and 8. The edge portion of the phase function 
controls the light rays that increase the light gathering and spatial 
resolution of the full aperture system but that without modification cause 
the largest amount of misfocus effects in traditional systems. It is these 
edge rays that should be modified most by the objective mask phase 
function because they control the variation of the MTFs and PSFs with 
misfocus. The actual modification made to these edge rays should position 
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them so that the sampled PSFs and MTFs are maximally insensitive to 
changes in misfocus. 

Notice that the MTFs from the Wavefront Coding Interference Contrast 
system of Figure 8 (upper and lower plots) essentially do not change with 
5 misfocus but also do not have the same shape as that of the in-focus MTF 

{\^ = 0) of the conventional Interference Contrast system. In the spatial 
domain, the Wavefront Coding Interference Contrast systems form Images 
with a specialized blur where the blur Is Insensitive to the amount of 
misfocus. The Image Processing function 326 ts used to remove this blur. 

1 0 The Image Processing function can be designed so that after processing the 

MTFs and PSFs of the combined Wavefront Coding Interference Contrast 
system, over a range of misfocus, closely match that of the in-focus 
Interference Contrast system. The Image Processing function can also 
produce an effective MTF that has more or less contrast than the in-focus 
1 5 Interference Contrast system, depending on the needs of the particular 

application. 

In essence, the image processing function restores the Wavefront Coding 
interference Contrast transfer functions to those expected from the 
conventional Interference Contrast system with no misfocus. Since all the 

2 0 Wavefront Coding MTFs are essentially identical, after image processing 326 

all MTFs (and hence all PSFs) will be nearly identical for each value of 
misfocus. 

More specifically, the image processing function, say F, implements a 
transformation on the blurred Wavefront Coding Interference Contrast 
25 system, say H^pQ so that after processing the system has an ideal 

response Hideai- Typically the Ideal response Is chosen as the in-focus 
response from the general Interference Contrast system. If Implemented as 
a linear filter, then F is (in the spatial frequency domain) equivalent to: 
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where w denotes a spatial frequency variable, (f the Ideal response is fixed 
then changing the Wavefront Coding Interference Contrast system H^^^ 
changes the image processing function F. The use of a different Wavefront 
Coding phase function can cause a change in the image processing function. 
In practice, It is common to be able to measure slight changes in the 
Wavefront Coding Interference Contrast system as a function of misfocus. 
In this case the image processing F is chosen as a best fit between the 
measured data and the desired system after processing. 

There are many linear and non-linear prior art techniques for removing 
known and unknown blur in images. Computationally effective techniques 
Include rectangularly separable or multi-rank linear filtering. Rectangularly 
separable linear filtering involves a two step process where the set of one- 
dimensional columns are filtered with a one-dimensional column filter and an 
intemiediate image Is fomied. Filtering the set of one-dimensional rows of 
this intennediate image with a one-dimensional row filter produces the final 
image. Multi-rank filtering is essentially the parallel combination of more than 
one rectangularly separable filtering operation. A rank N, digital filter kernel 
can be implemented with rectangulariy separable filtering by using N 
rectangularly separable filters in parallel. 

The form of the processing (rectangularly separable, multi-rank, 2D kernel, 
etc.) is matched to that of the Wavefront Coding element. Rectangularly 
separable filtering requires a rectangulariy separable Wavefront Coding 
element. The element described in figure 6 is rectangulariy separable. 

Figure 9 contains real worid images of human cervical cells made with a 
conventional Interference Contrast system and a Wavefront Coded 
Interference Contrast System. The image on the left of Figure 9 was made 
with a conventional 40X, NA=1.3 Interference Contrast system similar to 
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that of Figure 1 . The image on the right of Figure 9 was made with a 
Wavefront Coding Interference Contrast system similar to that of Figure 3. 
The Wavefront Coding Element 324 was a rectangularly separable cubic 
phase element. Rectangularly separable digital filtering was used for image 
processing 326. 

Notice the phase shading visible in the conventional image. This phase 
shading results in a 3D-like appearance of the object. This is a 
characteristic of Interference Contrast imaging. Notice also that many 
parts of the Interference Contrast images are blurred due to misfocus 
effects. The bottom part of the left image, for example, is particularly 
badly blurred by misfocus. The Wavefront Coded Interference Contrast 
image is also seen to have similar phase shading and 3D-like appearance as 
the conventional image. The depth of field visible In the image is much larger 
in the Wavefront Coded image than in the conventional image. Many parts 
of the cells that could not be resolved in the conventional image are clearly 
visible in the Wavefront Coding image. Thus, the Wavefront Coding 
Interference Contrast image produces both the characteristic Interference 
Contrast phase object imaging characteristics and a large depth of field. 

As shown in Figures 6 through 9, the Wavefront Coding Interference 
Contrast imaging system removes the effects of misfocus on the final 
Images. The Wavefront Coding Interference Contrast system will control the 
misfocus effects independent of the source of the misfocus. When increasing 
the depth of field, as shown in Figure 9, the misfocus effects are produced 
from the object or parts of the object not being in the best focus position 
relative to the imaging optics. Misfocus effects can also be produced by 
non-ideal optics, temperature changes/ mechanical positioning errors, and 
similar causes that lead to optical aberrations. Controlling misfocus effects 
besides those related to object positioning allows inexpensive systems to be 
produced that image with a high quality. For example, if the objective lens 
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312 of Figure 3 has a noticeable amount of chromatic aberration then 
misfocus effects will be produced as a function of illumination wavelength. 
The Wavefront Coding Interference Contrast system can control the 
chromatic aberration misfocus effects while also extending the depth of 
field. Other optical aberrations that can similarly be controlled include 
petzval curvature, astigmatism, spherical aberration, temperature related 
misfocus, and fabrication or alignment related misfocus. Any aberration 
that can theoretically or practically be improved in prior art systems by 
changing the focus position can also be conrected and Improved in 
Wavefront Coding Interference Contrast systems. 

What is claimed is: 
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CLAIMS 

1 . Apparatus (300) for increasing depth of field and controlling focus 
related aberrations in an Interference Contrast Imaging system having an 
Illumination source (302), optical elements (304, 306) for dividing light 
polarizations, and illumination optics (308) placed before a Phase Object 
(310) to be Imaged, and elements (314, 316) for recombining light 
polarizations and objective optics (312, 318) after the Phase Object to 
fonn an image (320) at a detector (322), the improvement comprising: 

an optical Wavefront Coding element (324) having an aperture and placed 
between the Phase Object and the detector, 

said coding element being constructed and an-anged to alter the optical 
transfer function of the Interference Contrast Imaging system in such a way 
that the altered optical transfer function is substantially insensitive to the 
distance between the Phase Object and the objective optics over a greater 
range of object distances than was provided by the unaltered optical 
transfer function, 

wherein the coding element affects the alteration to the optical transfer 
function substantially by affecting the phase of light transmitted by the 
element; and 

a post processing element (326) for processing the image captured by the 
detector by reversing the alteration of the optical transfer function 
accomplished by the coding element. 

2. The apparatus of claim 1 wherein the detector is a charge coupled 
device (CCD). 



3. The apparatus of Claim 1, wherein the phase of light transmitted by 
the coding element is relatively flat near the center of the aperture with 
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increasing and decreasing piiase near respective ends of the aperture. 

4. The apparatus of claim 1, wherein the phase of light transmitted by 
the coding element substantially follows a cubic function. 

5.. The apparatus of claim 4, wherein the phase of light transmitted by 
the coding element substantially follows a function of the form: 

Phase (x,y) = 12 [x3+ y^] 

where Ixl < 1 . lyl < 1 . 

6. The apparatus of claim 1, wherein the phase of light transmitted by 
the coding element substantially follows a rectangularly separable 
sum of powers function of the form: 

phase{x.y) = [ 2 a,- sign(x) Ixl ^+ C;Slgn(y) lyl ^/] 

where the sum is over the index /, 

sign{x) = -1 for x < 0, sign(x)=+l for x ^ 0. 

7. The apparatus of claim 1, wherein the phase of light transmitted by 
the coding element substantially follows a non-separable function of the 
form: 

phase(r,e) = [2 r^i cos( b, 9 + (j>j ) ] 
where the sum is again over the index /. 



8. The apparatus of claim 1, wherein the phase of light transmitted by 
the coding element substantially follows a function of the fonm: 
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Phase profile (x,y) = 

7[ sign(x) lx|3 + slgn(y) ly|3 ] + 7[ slgn(x) Ixl^-e + slgn(y) lyl^ ^ ] 
where Ixl ^1, lyl ^ 1 

9. The apparatus of claim 1, wherein the coding element further 
comprises a lens element for focussing the light. 

10. The apparatus of claim 1, wherein the coding element is integrally 
fomned with the Illumination optics. 

11. The .apparatus of claim 1, wherein the coding element comprises one 
of: 

an optical material having varying thickness; 

an optical material having varying index of refraction; 

spatial light modulators; or 

micro-mechanical mirrors. 
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12. The method for increasing depth of field and controlling focus related 
aberrations in a conventional Interference Contrast Imaging system having 
an illumination source, optical elements for dividing light polarizations, and 
illumination optics placed before a Phase Object to.be imaged, and elements 
for recombining light polarizations and objective optics after the Phase 
Object to fonm an image at a detector, the method comprising the steps 
of: 

between the Phase Object and the detector, modifying the wavefront of 
transmitted light; 

the wavefront modification step selected to alter the optical transfer 
function of the Interference Contrast Imaging system in such a way that the 
altered optical transfer function is substantially insensitive to the distance 
between the Phase Object and the objective optics over a greater range of 
object distances than was provided by the unaitered optical transfer 
function; and 

post processing the image captured by the detector by reversing the 
alteration of the optical transfer function accomplished by the element. 

13. The method of Claim 12, wherein the modifying step modifies the 
phase of light transmitted according to a profile which is relatively flat near 
the center of the aperi:ure with increasing and decreasing phase near 
respective ends of the aperture. 

14. The method of claim 12, wherein the phase of light transmitted by the 
element substantially follows a cubic function. 

15. The method of claim 12, wherein the phase of light transmitted by the 
coding element substantially follows a rectangularly separable sum of 
powers function. 
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16. The apparatus of claim 12, wherein the phase of light transmitted by 
the coding element substantially follows a non-separable function. 
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